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^ ■ Abstract 

,—1 ■ We consider quadrangulations with a boundary and derive exphcit expressions 

■ for the generating functions of these maps with either a marked vertex at a pre- 
I . scribed distance from the boundary, or two boundary vertices at a prescribed 

'. mutual distance in the map. For large maps, this yields explicit formulas for the 

2 I bulk-boundary and boundary-boundary correlators in the various encountered 

, scaling regimes: a small boundary, a dense boundary and a critical bound- 
ary regime. The critical boundary regime is characterized by a one-parameter 

^ , family of scaling functions interpolating between the Brownian map and the 

CN ! Brownian Continuum Random Tree. We discuss the cases of both generic and 

^ I self-avoiding boundaries, which are shown to share the same universal scaling 

tJ- ■ limit. We finally address the question of the bulk-loop distance statistics in the 

\^ ■ context of planar quadrangulations equipped with a self-avoiding loop. Here 

O , again, a new family of scaling functions describing critical loops is discovered. 

o : 

> : 

V~j • 1. Introduction 

><: 

■ Understanding the properties of random maps is a fundamental question for both 
the mathematical and the physical community. In mathematics, maps raise at a discrete 
level beautiful and rather involved combinatorial problems while, at a continuous level, 
they give rise to new probabilistic objects, like the Brownian map, whose construction 
is still under investigation. In physics, maps are used as discretizations for fluctuating 
surfaces in various domains ranging from low energy physics, for instance in the context 
of fluid membrane statistics, to high energy physics in the fields of string theory or of 
two-dimensional quantum gravity 0. 

The first incursion into these problems dealt mainly with global properties of ran- 
dom maps. At the combinatorial level, this amounted to a precise enumeration of various 
families of maps by several methods developed by mathematicians or physicists. These 
include in particular the original approach through recursive decomposition, developed 
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by Tutte in a series of papers 0, and the approach through random matrix integrals, 
which provide a systematic and powerful machinery for enumeration [0,3 . With this 
latter technique, the study was extended to maps with possible extra statistical degrees 
of freedom such as spins or particles (see for a review). At a continuous level, many 
results such as exponents characterizing global properties of the maps were obtained 
heuristically via the so-called Liouville model 0. Beside global properties, one then 
addressed the more refined question of the actual dependence of correlations on the 
distance along the map. A first expression was obtained in Refs. and (see also 
1^) for the so-called two-point function, which gives the law for the distance between 
two points on the map and more generally for the "loop-loop propagator", measuring 
the distance between two boundary loops at the extremities of a cylindrical map. These 
results were obtained in the context of triangulations via Tutte's recursive decomposi- 
tion approach, but at the price of heuristic arguments which, although non-rigorous at 
the discrete level, led eventually to the correct continuous correlators. Apart from this 
result, little remained known for quite a while on the statistics of distances in maps or, 
equivalent ly, on the metric structure of the Brownian map, probably because neither 
the recursive decomposition nor the matrix integral approach, nor even the Liouville 
model, are well-suited to address questions on the distance. 

Fortunately, a completely new enumeration technique was then discovered, where 
the distance plays a central role. It uses bijections to code the maps by much simpler 
objects such as the so-called well-labeled trees, where the labels precisely retain some 
of the distances in the map. This bijective approach was initiated by Schaeffer for 
quadrangulations (maps with tetravalent faces only) f^, and later extended to maps 



with arbitrary prescribed face valences [10|, Eulerian maps and maps with particles 



or spins As far as distance statistics is concerned, a first application dealt again 
with the two-point function, giving the law for the distance between two vertices on 



the map. This law was derived exactly at the discrete level in Ref. |]T2|] in the case of 
quadrangulations of fixed area, and was shown to converge to a continuous universal 
scaling function, giving a rigorous proof of the expression of Refs. and . A related 



quantity, the radius, was discussed in Refs. [13,14|. The universal continuous two-point 



function is an intrinsic characteristic of the Brownian map, a more fundamental object 
toward which many families of random planar maps (falling in the universality class 
of the so-called 2D pure gravity, such as maps with arbitrary bounded face degrees or 
maps coupled to non-critical statistical models) are expected to converge in the scaling 
limit where the area of the maps is large and scales as the fourth power of the distance 
1^,|16[. Transposing the Schaeffer bijection or its extensions at the continuous level 



allowed to construct this Brownian map as a random metric space, which was shown in 
particular to be homeomorphic to the two-dimensional sphere [|T7|,[T^] . Other properties 



of the Brownian map could be derived by first understanding their discrete counterparts 
and then taking a continuous scaling limit. For instance, the statistics of geodesies (i.e. 



paths of shortest length) between two points was considered in Ref. and it was 



shown that for typical points, all geodesies coalesce into a unique macroscopic geodesic 



path in the scaling limit [PC)| , |21[| . The so-called three-point function, which measures the 
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joint probability distribution for the pairwise distances between three uniformly chosen 
random vertices, was computed exactly in Ref. p2| in the case of quadrangulations 
of fixed area, and its universal scaling limit was analyzed. A remarkable property 



of confluence was discovered by Le Gall |21], stating that the three geodesic paths 



joining three typical points on the Brownian map merge by pairs before reaching their 
endpoints. Again, a complete characterization of the geometry of a geodesic triangle 
could be obtained from an exact solution at the discrete level [p3 |. 



So far, all the above rigorous results on the distance statistics obtained via the 
bijective approach dealt with closed planar maps, i.e with the topology of the sphere. In 
this paper, we extend these results to the more general context of maps with a boundary. 
More precisely, we are interested in the bulk-boundary correlator, which gives the law for 
the distance to the boundary of a vertex drawn at random in the bulk of the map, and in 
the boundary-boundary correlator, measuring the distance between two vertices on the 
boundary. Note that, in some sense, the bulk-boundary correlator lies half-way between 
the two-point function of planar maps and the loop-loop propagator of cylindrical maps. 
Our main results are explicit exact expressions for these correlators, already at a discrete 
level, in the particular case of quadrangulations. For quadrangulations of large area n, 
several regimes are obtained according to whether the boundary is small (with a length 
which remains finite), dense (with a length of order n) or critical (with a length of order 
7^1/2) We derive from our discrete results explicit expressions for the scaling limit of 
the bulk-boundary and boundary-boundary correlators in all these regimes. Most, but 
not all, of these scaling limits are universal and we recover in particular some results of 
Refs. and , here via a rigorous discrete enumeration. Maps with a critical boundary 
are characterized by a one-parameter family of universal scaling functions, corresponding 
to a new probabilistic limiting object interpolating between the Brownian map and the 
Brownian Continuum Random Tree. We first derive our expressions in the context 
of quadrangulations with generic boundaries, i.e boundaries which may have "pinch- 
points" separating the map into several components. Our results for the bulk-boundary 
correlator are then extended, both at the discrete level and in the various continuous 
regimes, to the case of self-avoiding boundaries where pinch-points on the boundary are 
forbidden. For large maps, generic and self-avoiding boundaries lead to the same scaling 
regimes. We finally address the problem of planar quadrangulations equipped with a 
self-avoiding loop. We give exact discrete and continuous expressions for the bulk-loop 
correlator, which is the law for the distance to the loop of a vertex drawn at random 
in the map. Here again, for large maps, a regime of critical self-avoiding loop is found, 
described by a new one-parameter family of scaling functions. 

The paper is organized as follows: Section 2 presents our main results, which are 
then proved in Sections 3 to 6. We give in Sect. 2.1 a precise definition of quadrangula- 
tions with a boundary and present in Sect. 2.2 a number of explicit discrete formulas for 
various generating functions encoding the bulk-boundary and boundary-boundary cor- 
relators. The case of both generic and self-avoiding boundaries are discussed. We then 
turn in Sect. 2.3 to the statistics of distances in large maps, for which several scaling 
regimes are found, depending on whether the boundary is small, dense, or critical. We 
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give in particular explicit expressions for various scaling functions characterizing these 
regimes. We finally discuss in Sect. 2.4 the bulk-loop correlator for quadrangulations 
with a self-avoiding loop, both at the discrete and continuous levels. Section 3 gives a 
precise derivation of our results for quadrangulations with a generic boundary at the 
discrete level. We present in Sect. 3.1 a bijection relating these quadrangulations to 
cyclic sequences of well-labeled trees. This property is used in Sect. 3.2 to derive ex- 
plicit formulas for associated generating functions. Physically, these correspond to the 
bulk-boundary and boundary-boundary correlators in a grand canonical ensemble with 
a boundary of fluctuating length, conjugate to a fugacity parameter z. We then deduce 
in Sect. 3.3 the corresponding fixed-length generating functions, corresponding to an 
ensemble of quadrangulations with a boundary of fixed length. Section 4 is devoted to 
the scaling limit of large maps with a generic boundary. We analyze in Sect. 4.1 the 
singularities of our discrete generating functions, which control the large map proper- 
ties. Three different scaling regimes are found: a small boundary regime for z < 1/8, 
discussed in details in Sect. 4.2, a dense boundary regime for z > 1/8, discussed in 
details in Sect. 4.3 and finally a "critical" regime for z ~ 1/8. This latter regime is best 
analyzed in Sect. 4.4 in the fixed length ensemble by considering large quadrangulations 
whose boundary has a length proportional to the square-root of their area. This gives 
rise to our one-parameter family of scaling functions. Section 5 deals with quadran- 
gulations with a self-avoiding boundary. We first show in Sect. 5.1 how to obtain an 
explicit formula for the corresponding bulk-boundary correlator from the expression of 
its generic boundary counterpart. This property is used in Sect. 5.2 to analyze the large 
map scaling limit, which is shown to present essentially the same three scaling regimes 
as above. Finally, Section 6 is devoted to quadrangulations with a self-avoiding loop. 
We show how to construct such quadrangulations by concatenating two quadrangula- 
tions with self-avoiding boundaries of the same length. This property is used to derive 
explicit expressions for the discrete and continuous bulk-loop correlators. In particular, 
a new family of scaling functions is found, describing maps with a critical self-avoiding 
loop. We end this paper by a few concluding remarks gathered in Section 7. 



2. Main results 

This section presents a panorama of our results, which are listed without derivation. 
We will then explain in the next sections how those can be obtained. 

2.1. Quadrangulations with a boundary: definitions 

Here and throughout the paper, we define a quadrangulation with a boundary as 
a planar map with a distinguished face such that all the other faces of the map have 
degree 4. Such maps are sometimes called "pseudo-quadrangulations" in the literature. 
We use the convention of representing the map in the plane with the point at infinity 
in the distinguished face (see Fig. 1), which we call the external face accordingly. Note 
that the degree of the external face is necessarily even. 
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Fig. 1: An example (a) of quadrangulation with a boundary of area n = 17 
and perimeter 2p = 34, which is the length of its contour (dashed green 
line). All faces are tetravalent, except for the external one which has 
degree 2p. Upon splitting the boundary at its separating vertices, we 
obtain (b) irreducible components (here six) which are either single edges 
or quadrangulations with a self-avoiding boundary. 

We call the boundary of the quadrangulation the set of edges and vertices incident 
to the external face. The actual sequence of edges followed when going around the 
external face counterclockwise in the plane is called the contour. The perimeter is the 
(even) length of this sequence, which is nothing but the degree of the external face. The 
bulk of the quadrangulation is the complement of the external face in the plane, and the 
number of inner faces it contains is called the area. 

As in the usual terminology, a pointed map is a map with a marked vertex, which 
is referred to as the origin, while a rooted map is a map with a marked oriented edge, 
which is referred to as the root edge. Here, we will make the useful convention that a 
rooted quadrangulation with a boundary has always its root edge on the boundary, with 
the external face incident to its right. We also consider in the following maps which are 
both pointed and rooted, or maps with are doubly rooted, in which cases all root edges 
are boundary edges with the appropriate orientation. 

Note that we consider here general maps which may possibly contain separating ver- 
tices or edges (a vertex or an edge is separating if its deletion disconnects the map). In 
particular, the boundary may contain such separating vertices or edges, which are those 
encountered several times along the contour. When the boundary has no separating ver- 
tices nor edges, it is said to be self-avoiding. Upon splitting at the separating vertices 
of the boundary, a quadrangulation with a generic boundary is naturally decomposed 
into several "irreducible" components, that are quadrangulations with a self-avoiding 
boimdary, or possibly single edges, arranged in a tree-like structure (see Fig. 1). This de- 
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composition will later allow us to study quadrangulations with a self-avoiding boundary- 
using results obtained for generic boundaries, which we present first. 



2.2. Discrete results 

We wish to enumerate quadrangulations with a (generic) boundary having pre- 
scribed area and perimeter. In the following, these two quantities will be usually de- 
noted respectively by n and 2p. As customary, the results are best expressed via a 
generating function, corresponding to a sum over all quadrangulations with a boundary, 
a given quadrangulation with area n and perimeter 2p having a contribution g'^z^. We 
will thus consider power series in two variables g and z. 

Arguably the simplest generating function is that for rooted quadrangulations with 



a boundary, already computed for instance in Ref. , and which may be written as 

Wo = W{1- gR^{W -1)) (2.1 
where R and W are the unique power series satisfying the algebraic equations 

R = l + 3gR^ 



W = 1 + zRW^. 



2.2) 



The particular form of this generating function yields, by two applications of the La- 
grange inversion formula, an explicit expression for the number of rooted quadrangula- 
tions with a boundary having area n and perimeter 2p: 

WoU^.. = 3-(2p)! (2n + p-l)! 

' pl{p - ly. n]{n + p + ly. ^ ^ 

where -Ig^^p denotes the extraction of the coefficient of g^z^ in the series. 

In this paper, we are interested in refined quantities involving the graph dis- 
tance. Our main results are exact expressions for the bulk-boundary correlator and 
the boundary-boundary correlator which are defined as follows. 

The bulk-boundary correlator is the generating function for pointed quadrangula- 
tions with a boundary, where the origin is at a prescribed distance, say d, from the 
boundary. It reads 

^ log (^) (2.4) 

where 

and X is the power series satisfying 

X = gR^{l + x + x'^). (2.6) 



6 



A few remarks are in order. The maps considered here might have non-trivial symme- 
tries, and the generating function Gd includes the corresponding usual inverse symme- 
try factor. Wd has itself a combinatorial interpretation, as the generating function for 
pointed rooted quadrangulations with a boundary, such that the origin is at distance 
smaller than or equal to d from the boundary, and such that the root edge starts from 
a (boundary) vertex closest to the origin. Such maps have no symmetries. Finally, the 
expression ( p.l|) is consistent with (|2.5| ) for d = 0, while Wd — > for (i ^ cxd, so that 
W is the generating function for pointed rooted quadrangulations with a boundary. 

The boundary-boundary correlator is the generating function for doubly- rooted 
quadrangulations with a boundary, such that the two root edges start from (bound- 
ary) vertices at a distance d from each other in the map. It reads 

Td = WiW - l)^/i (^-^ - 2{W - 1) + fd+iiW - If^ . (2.7) 

No symmetry factors are involved. 

Let us now consider quadrangulations with a self-avoiding boundary. A classical 
combinatorial argument shows that the generating function Wq for rooted quadrangu- 
lations with a self-avoiding boundary is related to Wq via 

Wo{9,Z) = Wo{9,z) with Z = zW^{g,z) (2.8) 

where we emphasize that Wq is a power series in two variables g and Z, Z being the 
variable conjugated to the half-perimeter of the self-avoiding boundary and g still being 
conjugated to the area. Algebra yields 

Wq^W{1- gR^iW -I)) (2.9) 

where R is still given by ( p.2| ) while W obeys 

W = 1 + ZR/{1- gR^iW -l)f. (2.10) 

By double Lagrange inversion we find the number of rooted quadrangulations with a 
self-avoiding boundary having area n and perimeter 2p 

W\ -o^-p (3p)! (2n + p-l)! 

Note that for self-avoiding boundaries we have the constraint 1 < p < n + 1. The 
bulk-boundary correlator, defined in the same way as for generic boundaries, reads 

Gd = log I ^ ^ — 1 for ci > 1, Go = Wq-1 (2.12) 

' Wd-i -Wo + l' ~ 
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where 



Wo 1-{W - l)fd 



Again the generating function Gd involves symmetry factors, unhke Wd which has the 
same combinatorial interpretation as Wd, now in the context of quadrangulations with 
a self-avoiding boundary. We have not been able to find a compact expression for the 
boundary-boundary correlator for self-avoiding boundaries. 

2.3. Distance statistics in quadrangulations with a boundary 

From these above exact expressions, we may now derive statistical information on 
distances in quadrangulations with a boundary. More precisely, we consider random 
quadrangulations with a boundary having prescribed area n and perimeter 2p, where 
each sample map appears with a probability proportional to its inverse symmetry fac- 
tor. We are particularly interested in the large n limit, for which we expect to find 
asymptotically the same results as with the uniform measure. Until further notice we 
consider generic (possibly non self-avoiding) boundaries. 

It proves convenient to consider first the fixed z ensemble where the area remains 
fixed equal to n but the perimeter fiuctuates, and each sample map with perimeter 
2p appears with probability proportional to (besides the symmetry factor). This 
model is well-defined (has a finite partition function for all n) for < z < 1/4. The 
above generating functions correspond to observables related to the distance. The bulk- 
boundary correlator encodes the distance between the boundary and a random vertex 
uniformly drawn in the bulk. The probability that this bulk-boundary distance be d is 
(G^d|g")/(log VF|g»), while the probability that it be less than or equal to d is nothing 
but [logWdlgn)/ {logW\gn). The boundary-boundary correlator encodes the distance 
between (the origins of) two edges uniformly chosen on the boundary. The probability 
that this boundary-boundary distance be d is {Td\gr^)/{2z-^Wo\gr^). 

For n — > oo, the model exhibits a phase transition at z = 1/8, which might be seen 
simply by analyzing Wq: 

- when 2; < 1/8 {subcritical regime), the perimeter remains finite as n — 00, 

- when ^ > 1/8 {supercritical regime), the perimeter is of order n and, up to Gaussian 
fluctuations of order n^/^, it concentrates around its mean value 2{p)n{z) with 

(P)n(^)~n^ (2.14) 

- when ^ ~ 1/8 {critical regime), the perimeter is of order n^/^. 

We will be especially interested in the critical regime, and for this case only we will 
perform the translation back to the fixed perimeter ensemble, considering quadrangula- 
tions with a boundary having fixed (large) area n and (large) perimeter 2p, keeping the 
renormalized half-perimeter P = p ■ n~^/^ finite. Let us now discuss the manifestations 
of the transition on the distance statistics, as seen by analyzing Gd and Td- The follow- 
ing table gives a qualitative summary of the asymptotic behaviors for the perimeter, the 
bulk-boundary distance and the boundary-boundary distance in the various regimes. 
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perimeter 


bulk-boundary 
distance 


boundary-boundary 
distance 


universality 
class 


z < 1/8 


finite 


^1/4 


finite 


Brownian map 


z > 1/8 


n 


finite 


^1/2 


Brownian CRT 


z ~ 1/8 


^1/2 


nV4 


nV4 


Brownian map with 
a boundary 



In the subcritical regime, the bulk-boundary distance is of order n^/^ and admits a 
continuous hmit law which does not depend on z (provided it is in the subcritical range 
]0, l/8[), and whose (cumulative) distribution function reads 



lim 

n — >oo 



log 



logW\ 



9 /"CXD 

di\ie-^":F{D--e) (2.15) 



where 



\ 



1 + 



\ 



sinh 



1/41) 



(2.16) 



We recognize the universal two-point function of pure 2D gravity ^ also obtained 
in the case of quadr angulations without a boundary [0. This result agrees with the 
physical intuition: upon rescaling distances by a factor n~^/^, the boundary reduces to a 
point which behaves no different from a typical point in a large random quadrangulation. 
Mathematically, the metric space obtained in the scaling limit is expected to be the 
Brownian map. In contrast, the boundary-boundary distance remains finite and admits 
a non-universal discrete limit law with a computable, albeit complicated, expression. 

In the supercritical regime, the bulk-boundary distance remains finite at large n 
and admits a non-universal discrete limit law, whose cumulative distribution function 
reads 



■Xd) 



, logVFdl^. 
lim 



where 



3^crit(^) — 



{1 - {x,,-Az)r+'){i - {x,rAz)r+') 

162 - 1 - v/3((8z)2 - 1) 



(2.17) 



(2.18) 



2(1 -4^) 

The boundary-boundary distance is of order n^/^ and admits a continuous limit law. 
Using the square root of the mean half-perimeter ( ^.141 ) as distance unit, we find the 
Rayleigh probability density function 



Pbound. (5) 



lim 

n — >cxD 



T, 



2z:fWo\ , 

az ^ I q' 



V {P)n{z) 



25 e 



(2.19) 
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which coincides with the two-point function for the Brownian Continuum Random Tree 
||25|| . The (Brownian) CRT is expected to be the hmiting metric space obtained when 
rescahng distances by a factor n~^/^. The physical interpretation is that, in the super- 
critical phase, the boundary becomes "dense" in the quadrangulation, and folds onto 
itself, creating a branched structure. 
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D 



Fig. 2: Plots of the cumulative distribution function $(-D, P) as a function 
of D, for P = 0.01, 0.1, 0.5, 1.0, 2.0 and 5.0 (thin lines from bottom to 
top) . We also indicate in thick red line the (integrated) two-point function 
$(D). 

In the critical regime, both the bulk-boundary and boundary-boundary distances 
are of order n^/^. Both admit continuous limit laws, which are best expressed in the 
(critical) fixed perimeter ensemble. On the one hand, the bulk-boundary distance cu- 
mulative distribution function reads 



$(D, P) = lim 

n — >oo, pr^Pn^/"^ 



log W\ nr,^/4. I I 



where 

n{p,p-ix) 



/(I^;^) = y|//^ coth 



= e^'/4 r di ^ e-«' n{D, P; -C') 

J-oo 1 

/•oo 2 



(2.20) 



(2.21) 



It is plotted in Fig 2 as a function of D, for P = 0.01, 0.1, 0.5, 1.0, 2.0 and 5.0. As 
we will show later, for P — we recover the two-point function $(-D) (shown in red 
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bound. 



. 6 



. 4 



0.2 



P = 




P=0 



Fig. 3: Plots of the probability density pbound.('^, P) as a function of the 
rescaled distance 5 = D/y/P for P = 0.5, 1.0, 1.5, 2.0, 3.0, 5.0, 10.0 
(thin lines from bottom to top). We also indicate the Rayleigh law of 
Eq. ( |2.19| ) (upper blue thick line) corresponding to the limit P — > oo and 
the non-trivial law of Eq. (2.25) (lower red thick line) corresponding to the 
hmit P ^ 0. 

on Fig 2), while for P — oo the rescaled bulk-boundary distance is of order P~^/^ and 
$(D, P) takes the simple scaling form: 



$(P',P) ~ tanh^ I /P 



(2.22) 



consistent with the supercritical law ( p.l7|) for z 1/8"*". On the other hand, the 
boundary-boundary distance probability density function reads 



PhoundXD,P) = lim 



(7"2P 1/4 

■ n ' 



3P4 



g-D^/P |(2P)3_3i}p) + (4i)2p_ 2p2)^^(i}^ P)+ 2DPV2(P', P)} 

(2.23) 



where 



ai(D,p) = ^ / ^ e-^'+^^^m-e 

VnP 1 



<T2{D,P) 



(2.24) 
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and f{D,fj,) is given by ( |2.21| ). It is also natural to measure the boundary-boundary 
distance in units of ^yp which amounts to introducing the variable 5 = D/^/P. The 
corresponding probability density pboundX'^' is plotted in Fig. 3 for P = 0.5, 1.0, 1.5, 
2.0, 3.0, 5.0, 10.0. When P — > oo we precisely recover the Rayleigh density ( p.l9|) , while 
for P ^ we have the particularly simple but non-trivial expression: 

Pbound.{S,pf^''^e-'\356 + 28S^ + 126'' + 36^) . (2.25) 

We emphasize that all these expressions are expected to be universal (up to a 
possible rescaling of the distance and perimeter) and are intrinsic to the metric space 
obtained in the scaling limit. More precisely, up to a change of the distance scale, we 
have a one-parameter family of random metric spaces indexed by the (renormalized) 
perimeter P. It might be called the Brownian map with a boundary. Note that P is not 
homogeneous to a distance but to its square, an indication that the fractal dimension 
of the boundary is two. The Brownian map with a boundary interpolates smoothly 
between the Brownian map, recovered for P = 0, and the Brownian Continuum Random 
Tree, recovered in the limit P — > oo. Note that when P 0, we observe a deviation 
from the Brownian map statistics for small distances of order \/P (corresponding to 
finite values of 5 above). 

Let us now briefly mention the results for quadrangulations with a self-avoiding 
boundary. Observe that the generating functions Wq and Gd ai'e related to random 
quadrangulations with a self-avoiding boundary in exactly the same way as Wq and Gd 
are related to random quadrangulations with a generic boundary. A parallel approach 
can thus be followed. Here Z denotes the activity per unit of half-perimeter. There 
is now a phase transition at Z = 2/9 which is expected to be in the same universality 
class as the above. From the exact expressions for Wq and Gd we can show that: 

- for Z < 2/9, the perimeter remains finite as the area n tends to infinity, while the 
bulk-boundary distance is of order n^/"* with the same limit law ( |2.15| ), 

- for Z > 2/9, the perimeter is of order n and concentrates around its non-universal 
mean value, while the bulk-boundary distance is finite and has a non-universal 
discrete limit law, 

- for Z ~ 2/9, the perimeter is of order n^/^, while the bulk-boundary distance is of 
order n^^^. 

In this latter case, we may as well consider the critical fixed perimeter ensemble, and 
compute the bulk-boundary distance cumulative distribution function. We find the 
same expression as in (|2.20| ) up to a factor 3 in the renormalized perimeter P (see 
details in Sect. 5.2). This is a first non-trivial check of the universality of our analytical 
expressions. 

2.4- Application to self- avoiding loops 

Another interesting application of our exact discrete results is that they allow us 
to study the statistics of distances in quadrangulations with a self-avoiding loop. More 
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precisely, a self-avoiding loop is a closed path made of consecutive edges of the quad- 
rangulation, which is simple, i.e visits any vertex at most once. We consider planar 
quadrangulations with a distinguished oriented self-avoiding loop (and no boundary: 
all faces have degree 4). The area is the total number of faces, while the loop length is 
necessarily even. Again we consider a statistical model where the area n is fixed, where 
the loop length 2p may either be fixed or be controlled by a weight y^, and where in all 
rigor we need to incorporate the inverse symmetry factor, irrelevant for ti ^ oo. This 
is a particular instance of the so-called ©(A/" = 0) model on a random lattice, slightly 
different from the ones studied with matrix model techniques where the loops would 
run on the dual map. 

The connection with quadrangulations with a boundary is easily seen. Upon cutting 
along the loop, a quadrangulation with a self-avoiding loop yields two quadrangulations 
with a self-avoiding boundary, constrained to have the same perimeter. The orientation 
of the loop allows to distinguish these two pieces as left and right, and we clearly have 
a bijection preserving the total area. We can therefore express a number of generat- 
ing functions for this problem in terms of the generating functions found above. For 
instance, the generating function for quadrangulations with a self-avoiding loop and a 
marked vertex on the loop reads 

^oig.y) = J2y' {Woi9.Z)\z.y (2.26) 
P>1 

where g is the weight per face while is the weight per edge of the loop. More 
generally, we may consider the generating function for quadrangulations with a self- 
avoiding loop and a marked vertex at distance d from the loop and lying on its right. 
It reads 

Td{g,y) = J2y''Wo{9,Z)\z.Gd{9,Z)\z. (2.27) 
p>i 

if the configurations are counted with their inverse symmetry factor. Constraining the 
marked vertex to be on the right of the loop ensures that both expressions are consistent 
for (i = 0, and by symmetry it causes no loss of generality. Statistically, encodes 
the bulk-loop distance, i.e the distance between the loop and a random vertex uniformly 
drawn in the bulk. In the following sections, we provide more explicit (yet slightly 
involved) expressions for Wo\zp and Gd\zp, easing the task of deducing the bulk- loop 
distance statistics for maps of large fixed size n. To sum up our results, we find a phase 
transition at y — 4/81. 

- For y < 4/81, the loop length remains finite as n ^ oo, and the bulk-loop distance 
is of order n^/^ with a distribution again characterized by the two-point function 
( p.l5|) . The physical interpretation is that the loop remains microscopic and is thus 
irrelevant in the scaling limit, where distances are rescaled by a factor n~^/^, and 
which is still described by the Brownian map. 

- For y > 4/81, the loop length is of order n and the bulk- loop distance is finite. The 
physical interpretation is that the loop becomes dense in the quadrangulation. We 
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however lack evidence that the scahng hmit (on a scale n^/^) is still described by 
the Brownian CRT, though this hypothesis is plausible. 
- For y ~ 4/81, the loop length is of order n^/^ and the bulk- loop distance is of order 




D 



Fig. 4: The cumulative distribution function $(1), P) as a function of 
the the bulk-loop distance for a self-avoiding loop of (rescaled) half-length 
P = 0.01, 0.1, 0.2, 0.5 and 1.0 (thin lines from bottom to top). We also 
indicated (thick red line) the limiting two-point function $(-D). 

Here again we are most interested in this critical case, and the results are best 
expressed in the ensemble where both the area n and the loop length 2p are prescribed, 
and jointly taken to be large keeping the ratio P = p ■ n~^^'^ finite. The scaling law for 
the bulk-loop distance cumulative distribution function reads 



Hd,p) 



lim 

n — >oo, pr^Pn^/"^ 



fc=0 k\gn.yV 



(2.28) 



18^^^ 



1 + 18P2 



r di ^ e-^" HiD, P; -e) 

J — oo 1 



where 

n{D,P-,fi)^3-^^^{mP^)ll + {3^-2f{D-,fi)) rfKe-^-2^(^'^)^2K 

(2.29) 

It is plotted in Fig 4 for P = 0.01, 0.1, 0.2, 0.5 and 1.0. When P 0, $(£>, P) ^{D) 
as expected while, when P oo, D is of order P~^/^ and we have the scaling form 



$(P',P) ~ tanh^ 




(2.30) 
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The scaling function P) is expected to be universal, and characteristic of a model 
of self-avoiding loop on a Brownian map. 



3. Quadrangulations with a boundary: combinatorics 

We now come to the derivation of the expressions given in Sect. 2.2 for the var- 
ious generating functions concerning quadrangulations with a generic boundary. Our 
approach is based on a bijection with simpler objects, namely sequences of well-labeled 
trees, as discussed just below. 

3.1. Bijection 

A quadrangulation with a boundary is a particular instance of a bipartite planar 
map. As such, it may be coded by a so-called well-labeled mobile, as explained in the 
section 2 of Ref. |Ty] (see Fig. 5). More precisely, the coding is for a pointed map (i.e a 
map with a chosen origin vertex). The associated mobile is a plane tree with alternating 
labeled and unlabeled vertices. The labeled vertices correspond to the original vertices 
of the map and they carry an integer label equal to the graph distance in the map from 
the corresponding vertex to the origin. The unlabeled vertices of the mobile correspond 
to the faces of the map and their degree is half the degree of the corresponding face in 
the map. Around each unlabeled vertex v, we have the following property (P): reading 
the sequence of labels of vertices adjacent to v clockwise around v, any label I is followed 
by a label larger than or equal to £ — 1 . Finally, the mobile has a minimum label equal 
to 1. 

For quadrangulations with a boundary of length 2p, all the unlabeled vertices of 
the associated mobile necessarily have degree two, except for that associated with the 
external face, which we call the external vertex and whose degree is p. Those bi-valent 
unlabeled vertices may be erased, giving rise to edges which connect the labeled vertices 
directly. The property (P) may then be rephrased into the property (P') that labels 
on adjacent labeled vertices differ by at most 1. The resulting object is therefore a 
collection of p well-labeled trees, i.e trees with labeled vertices satisfying (P'), attached 
to the external vertex by their root vertices, whose clockwise sequence of labels satisfies 
(P) around the external vertex (see Fig. 6- (a)). An equivalent but more convenient 
coding of the sequence of these p root labels around the external vertex is via a cyclic 
sequence of 2p non-negative integers such that consecutive integers differ by ±1 (see 
Fig. 6-(b)). In this coding, the root labels simply correspond to those integers which 
are followed immediately by a smaller integer in the cyclic sequence, and the property 
(P) is automatically satisfied. Moreover, the cyclic sequence of integers corresponds 
precisely to the distance to the origin of the successive boundary vertices along the 
contour (see Fig. 6-(b)). 

To summarize, we have a bijection between, on the one hand, pointed quadrangu- 
lations with a boundary and on the other hand cyclic sequences of non-negative integers 
such that consecutive integers differ but ±1, with a well labeled tree with root label i 
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Fig. 5: An example (a) of pointed quadrangulation witli a boundary. We 
have labeled each vertex by its graph distance from the origin (vertex with 
label 0) . Adding (b) in each face an unlabeled vertex and connecting it to 
those labeled vertices followed by a smaller label clockwise within the face, 
we end up (c) with a particular labeled mobile whose unlabeled vertices all 
have degree 2, except for that associated with the external face, which has 
a degree equal to half the perimeter of the quadrangulation. The labels 
around this vertex satisfy the property (P) of the text. 

attached to each descending step £ ^ £—1 of the cyclic sequence, and with the require- 
ment that the global minimum label is 1. Under this bijection, we have the following 
correspondences : 

- The total number of edges for all the well-labeled trees is equal to the area n 
(number of inner faces) of the quadrangulation. 

- The label of any vertex is equal to the distance to the origin of the corresponding 
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Fig. 6: The well-labeled mobile of Fig. 5-(c), where we erased the bi- 
valent unlabeled vertices (a), may alternatively be viewed (b) as a set of 
well-labeled trees satisfying the property (P') of the text, attached at each 
descending step of a (counterclockwise-oriented) cyclic sequence of integers 
(green outer circle) reproducing precisely the distance from the origin of the 
successive vertices along the contour of the quadr angulation (as apparent 
in Fig. 5- (a)). 

vertex on the map. 

- The (even) length of the cyclic sequence is equal to the length 2p of the boundary. 

- The successive integers in the cyclic sequence are equal to the distance to the origin 
of the successive boundary vertices along the contour. In particular, the smallest 
integer d in this cyclic sequence is the smallest distance found between the origin 
and a vertex of the boundary, i.e the distance from the origin to the boundary. 

We may transform the cyclic sequence of non-negative integers into a Dyck path 
(0 = ■ ■ ■ 5 ^2p-i5 ^2p = 0) by reading the sequence from one of its minima, and 

subtracting d from all the integers. To each descending step ii ii — 1 is now attached 
a well-labeled tree with root label ii + d. Note that the cyclic sequence may have sev- 
eral minima, and choosing a particular minimum amounts to marking a boundary edge 
d ^ d + 1 with, say, the external face to the right. 

3.2. Basic generating functions 

We now wish to compute the generating function for pointed quadrangulations 
with a boundary, where a map with area n and perimeter 2p comes with a weight g'^z^. 
Under the bijection, this weight simply amounts to a weight g per edge of the well- 
labeled trees and a weight z per descending step of the cyclic sequence of integers, or 
equivalently of the Dyck path. 

A first ingredient is the generating function Ri for rooted well-labeled trees with a 
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root label i > 1 and with the condition that all the labels on the tree are larger than or 
equal to 1. This generating function was computed in Ref. fl^ and reads 

'I [^ + 3]. 



Ri = R 



[i+ll [i + 2l 



where we use the notation 



1 — X 
1 — X 



and where the quantities R and x are solutions of 

R = l + 3gR'^ , X = gR^{l + x + x^) 

namely 



R 



X 



1- y/T^n^ 

6g 

1 - 24(7 - yi - 12(7 + v^V72^2 + Qg + ^l-12g-l 
2(6(7 +v/T^n^-l) 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



min-^'('u) > 1 




Fig. 7: Opening a cyclic sequence of non-negative integers at one of its 
minima d results into a Dyck path (0 — £o, £i, ■ ■ ■ , l2p-i^ ^2p = 0) by sub- 
tracting d from all integers in the sequence. To each descending step 
£ £ — 1 of the Dyck path is attached a well-labeled tree with root label 
I + d, resulting in a weight zRd+i for the descending step if we demand 
that all the original labels be larger than 1. 

We can now express the generating function Wd for Dyck paths with a well-labeled 
tree with root label £i + d attached at each descending step £i — > £j — 1, and with a 
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global minimum label larger than or equal to 1. This generating function reads 

Wd = Y. E n ^^^^+^ (3.5) 

Dyck paths of length 2p descending steps 

- (0=lo,l'i^,...,l2p=0) 

with a conventional weight 1 for the trivial Dyck path of length 0. We have the recursion 
relation 

Wd = l + zRa+iWaWa+i , (3.6) 

as obtained by decomposing any non-trivial Dyck path into its first ascending step 
^ 1, a path from 1 to 1 with all intermediate heights larger than or equal to 1 (weight 
Wd+i as obtained by a simple shift of labels), its first descending step 1^0 (weight 
zRd+i) and a final Dyck path (weight Wa)- 

We may look for a solution of this equation in the form 

^^ = ^i7Tlr^' (3.7) 



with W and Va to be determined. Substituting this particular form in Eq.( |3.6|) leads to 
the equation 

W[d + 2]^ Vd+i = [d + 3], Vd + zRW"" [d + 1], Vd+2 . (3.8) 
This relation is satisfied (for arbitrary R and x) upon taking 

Vd = l + Xx'^ (3.9) 

provided we choose W and A such that 



W = 1 + zRW'^ 



Wi-x"^ + Ax) = {-x^ + A) + zRW^{-x + Xx^) 

namely 



(3.10) 



2zR ' l-x{W-l) ^ ' 



Plugging this last expression back in (|3.9| ) and (^j) ' with the desired solution 

[d + 2]^ [rf + 3],-x(M^-l)[d+lL 



= W 



[d + K [d + 2\^-x{W-l)[dl 
l-{W-l)fd+i 



(3.12) 



l-{W-l)fd 
where W is explicitly given in terms of g and z as 



(1 - .fr^)z ^ ' 
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while we introduce the compact notation 



with l']^ defined above. In particular, we have fo = and /i = gR^. Note that Wd 
depends on the variable z only through the quantity W. In particular, for d = 0, we 
find 

Wq = W{1- gR^{W -1)) (3.15) 
which agrees with the known generating function for rooted pseudo-quadrangulations 

To go from Dyck paths back to cyclic sequences of integers, we simply have to 
identify the Dyck paths differing only by the choice of an instance of the smallest 
integer in the sequence. We note that Eq. ( |3.6|) may be alternatively written as 



Wd = J2i^Rd+iWd+i)'' (3.16) 

fc>0 

where (zRd+iWd+i)'' is nothing but the generating function for Dyck paths with exactly 
k returns to 0. Cyclic sequences are then enumerated by 

J2 UzRd+iWd+i)'' = logWd . (3.17) 

k>i 

where we use the convention of counting configurations with an inverse symmetry factor. 
A given configuration made of a cyclic sequence with its attached trees may only have 
a cyclic symmetry group Z^, (with m a divisor of the number of minima k), and is then 
counted with a weight 1/m. 

So far, in the configurations counted by logWd, we imposed only that the global 
minimum label be larger than or equal to 1. To impose that this minimum label be 
exactly 1, as required by the bijection, we must suppress those configurations with a 
minimum larger than or equal to 2, whose generating function is obtained from the 
previous one by a simple shift of all labels by —1, i.e is given by log Wd-i- In particular, 
we deduce that the generating function Gd for quadrangulations with a boundary and 
a marked vertex at distance d from the boundary is simply given by 

G.^.„g(^). (3.18) 

Again, configurations with an m-fold symmetry around the marked vertex are counted 
with a factor 1/m. 

As for the original quantity Wd, it is the generating function for quadrangulations 
with a boundary, with a marked vertex at a distance less than or equal to d from the 
boundary (the origin) and with a marked "closest edge" to this origin, i.e a boundary 
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edge incident to a vertex at distance d from the origin and oriented counterclockwise 
around the bulk of the quadrangulation. Note that such pointed rooted maps cannot 
have any non-trivial symmetry. Finally, we may also interpret log Wd as the generating 
function for pointed quadr angulations with a boundary whose origin is at a distance 
less than or equal to d from the boundary. 

As discussed in Section 2, we are also interested in the generating function Td for 
quadrangulations with a boundary having two marked (and distinguished) boundary 
edges oriented, say counterclockwise around the bulk of the quadrangulation, such that 
the origins of these marked edges are at a mutual distance d on the map. Taking the 
origin of the first marked edge as the origin of the map, we get under the bijection a cyclic 
sequence of integers with minimal value 0. The first marked edge defines a first step 
1 at which we may start reading the cyclic sequence, leading to a Dyck path with a 
well-labeled tree with root label ii attached to each descending step ii — Marking 
the second boundary edge amounts to choosing a step d ^ d±l in the sequence, i.e to 
choosing a point of height d in the Dyck path (for d = 0, this point must be different 
from the last one). Upon decomposing the Dyck path into a first ascending part from 
to the marked point d and a second descending step from (i to 0, we get the expression 

Td = WoWi---Wdx{WdzRd){Wd-izRd-i)---{WizRi)Wo 

( o. 19) 

= {WoWi---Wdfz''{RiR2---Rd) , 

valid for d > 0. For d = 0, the actual decomposition yields Tq = Wq — Wq. From now 
on, we will assume d > when referring to Td. Upon substituting the expression ( |3.12|) 
for Wd, we get 



,2 



[3]. [^+l].[^^ + 3]. ^3_2o) 

= W^iW - ifh (^-^ - 2{W - 1) + fd+i{W - If^ . 

3. 3. Fixed length generating functions 

As we already noticed, the quantities Wd and Td depend on z only via the quantity 
W as given by Eq. ( |3.1C1| ), which we may rewrite as 

(3.21) 



i?(l + «;)2 
upon introducing the notation 

w = W -I . (3.22) 
In this respect, Wd and Td are so-called Lagrangean generating functions, i.e for which 



we can apply the Lagrange inversion theorem [B^. This means that we may extract an 



explicit expression for the z'p term of these generating functions. More precisely, the z'p 
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term, expressed as a contour integral in the variable z around may be transformed by 
the change of variable z ^ w into a contour integral around of the variable w, namely 



dz 1 



RV ^ ^_L^(1_^)(1+^)2P-1 



2i7r zP+i I J / 2i7r w^+i 

Upon expanding ( p.l2|) in w = 1^ — 1, we get the expression 

/ 

t^d = (! + «') 

V 

Taking the contour integral (|3.23|) of this quantity, we immediately get 



(3.23) 



(3.24) 



2p\ f 2p 



P J \P~ 
{2p)\ 



-A -{fd+i-fd)^ (( 
' k>i 



2p 



2p 



p — kj \p — l — k 



k-1 



ifd+l-fd) ^ 



{2py. 



p\{p+i)\ ^^{p-k)\{p+k+i)\ 

By a similar calculation, we easily obtain 



{2k + l){fa) 



k-i 



(3.25) 



(logl^d)l 



zP 



for p > 1. For d = 0, we get in particular 

with R given by (^^- We immediately deduce that 

3''(2p)! {2n + p-l)\ 



(3.27) 



W, 



0\g"zP 



p\{p — ly. n\{n + p + ly. 
for p > 1. On the other hand, for d oo, we have Wd W, with 

(2p)! 



(3.28) 



so that 



3''{2py {2n + p-iy 

~ {p- iy{p + iy ni{n+py 
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(3.29) 



(3.30) 



for p > 1. The quantity Wo|g"zP is the number of quadrangulations with area n and 
boundary of length 2p, where we have marked one of the 2p boundary edges on the 
contour. The quantity Wlg-n^p is the number of the same quadrangulations with addi- 
tional markings of one of their n + p + 1 vertices (at some arbitrary distance from the 
boundary) as well as of a particular boundary edge closest to this vertex. The quantity 



Wlg-nzP 2p 2p 



WolqTz^p n + p + 1 p + 1 



(3.31) 



measures therefore the average number of boundary edges closest to a uniformly chosen 
random vertex, for the ensemble of quadrangulations with area n and perimeter 2p. 
Note that, surprisingly, this average number is independent of n (and its value thus 
matches that obtained for ti = 0, where quadrangulations with a boundary reduce to 
plane trees with p edges). 

Finally, we may write ( |3.20|) as 

T, = (1 + «;)2/i - 2«;"+^ + (3.32) 



so that we get 
Td\zp =R^fi 



1 ff2p + 1\ f 2p + l 



fd+i \\ p-d J \p-l-d 
p — l—dj \p — 2—dJJ \\p—2—dJ \p — 6—d 



2 op. [ (c^+l) (2p + l)! 

{^p-d)\{p+d+2)\ 



(3.33) 



'{p-d-l)\{p+d+Z)\^ '■"'^\p-d-2)\{p+d+A)\ 

By a slight refinement of this calculation, we can compute the generating func- 
tion Trf(s,s') for quadrangulations with a boundary of length 2p = s s' with two 
marked boundary edges whose origins are separated by s steps counterclockwise along 
the contour and s' steps clockwise, and are at a mutual distance d on the map. We find 

Td{s- s') = R'+O, I -r^C(s, d)C{s', d) - C{s, d)C{s\ d + 2) 

^Jd+i (3.34) 

- C(s, d + 2)C(s', d) + fd+iC{s, d + 2)C{s', + 2) } 

where 

• (3.35) 

This results holds when d, s and s' have the same parity, while Td(s, s') vanishes other- 



wise. 
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4. Quadrangulations with a boundary: asymptotics 



4.1. Critical lines 

In order to describe ttie statistics of distances in large quadrangulations, we have to 
analyze the singular behavior of the various generating functions above. More precisely, 
the asymptotics for a large area n is encoded in the singularity reached at the radius 
of convergence in g of these generating functions. It is simpler to work first with a 
fixed value of z, corresponding to the fixed z ensemble mentioned in Section 2. A first 
singularity is associated with the singular behavior of R and x, as given by Eqs. ( |3.3|) 



or (|3.4|) , when g approaches the critical value 

ail = ^ (4.1) 

irrespectively of the value of z. This is the dominant singularity for z < 1/8 while, 
for z > 1/8, another singularity comes from the singular behavior of W, as given by 
Eqs. (13.101 ) or (|3.11| ), when zR approaches the value 1/4, which defines the critical line 

^7£iW = ^4;2(l-4z) (z>l/8). (4.2) 

(2) 

Since (/criu-^) — 1/12, this second singularity is dominant (i.e determines the radius of 
convergence) whenever z > 1/8. 

The radius of convergence therefore changes determination at 

^crit = ^ • (4.3) 

As we shall see below, the generating functions have very different scaling behaviors 
when z is smaller or larger than ^crit- As discussed in Section 2, this is the manifestation 
of a drastic change in the geometry of large quadrangulations with a boundary at this 
critical value. 

4-2. Scaling limit: the z < 1/8 regime 

Let us first discuss the case ^ < 1 /8 for which the generating functions have radius 
of convergence g^l^^ = 1/12. We may analyze the associated singularity by setting 

9=^{l-^^e) (4.4) 

with e ^ 0. A sensible scaling limit for Wd is obtained by considering large distances 
of the form 

d = De-^/\ (4.5) 
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with D finite. We then have the foUowing small e expansions 

fa = l- 2f{D; ^) e^^ + {Af\D- ^) - 2^) e''^ + ■ 
fd+i = 1 - 2f{D- ^) e'/^ + (6/2(D; ^) - 6^;^) e'/^ + ■ 



where we define 

We also have the expansion 

W = A{z) - zA'{z)y/iit^'^ + 

where we define 



A{z) 



1 - VI -8^ 



Plugging these expressions in ( p.l2|) , we obtain the expansion 

Wd = A{z) - T{D; ^) z J^{z) + • • 
where we introduce the notation 



(4.6) 

(4.7) 

(4.8) 
(4.9) 

(4.10) 



sinh^ (^^1 //i/4Dj 



(4.11) 



We may alternatively derive ( [4.1C1| ) by using the explicit form ( |3.25| ) of VFdIzP and ex- 
panding it at small e. Using the expansion 



(4.12) 



we obtain that 



(2p)! 



p\{p+l) 



.1/2. 



{2p)\ 



(p-l)!(p+l) 



-{2fiD;^,)-3^)e'/'Yl 
{2pV 



{2p)l 



2P- 



pl{p+ 



j^J^ {p-k)\{p+k + l)\ 
l-{l-p^(D;^)ei/2 + ---} . 



(2A; + 1) } + 



(4.13) 



Upon summing over p with a weight z^, this reproduces precisely the expression (|4.10| ). 
Note that the "kernel" J-'{D; fx) occurring in (|4.13|) is independent of p. For p = 1, quad- 
r angulations with a boundary of length 2 are equivalent, upon closing the boundary, to 
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rooted quadrangulations, i.e quadrangulations with a marked edge. The kernel J^{D; fx) 
is therefore the same as that encountered in Ref. [jl2| when deriving the continuous 
two-point function of planar quadrangulations from the generating function of rooted 
quadrangulations. 

We may repeat this analysis for the quantity \og{Wd)- We have the expansion 

log Wd = log Aiz) - J^{V; /.) + • • • (4.14) 

or the equivalent expansion 

(logWa)U.=2^^^^^{l-pnD;f,)e'/' + ...] . (4.15) 

From the above singularity analysis, we may deduce the bulk-boundary distance 
statistics when 2; < 1/8 in the ensemble of quadrangulations (with a boundary) with a 
fixed area n, in the limit where n 00. Indeed, we may extract the contribution to 
Wd or log Wd of these quadrangulations by a contour integral around in the variable 
g. At large n, this translates into an integral over a real variable ^ upon setting (see 
Ref. Ijl^ for a more detailed discussion) 



9=^1 + ^). (4.16) 



We may indeed write at large n 



Setting 

d = Dn^/'^ (4.18) 

with D finite, we may use the expansion ( [4.10| ) above with e = 1/n and fi = — and 
deduce that, at large n 

Wd\gr^ JH^^i^) / d^i^, e-^ T{D--e) (4.19) 

since the first (regular) term A[z) in the expansion ( [4.1UD leads to a vanishing integral 
in ^ by parity. In particular, for L> — 00, we have J-'{D,—^'^) — > — i^ and we get 

12"- 

which may alternatively be obtained directly from the general expression (|3.30|) . 
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Similarly, we obtain from ( |4.15|) that 

rn-^/^ A{z) 



\ogW\ 



IT" zA'{z) 



(4.21) 



2v^?l3/2 A{z) ' 
The ratio of these quantities tends at large n to a finite quantity 

^{D) = — d^i^ e-^V(I); -e) (4.22) 

which is the (cumulative) distribution function for D giving, in the ensemble of pointed 
quadrangulations with a boundary, the probability that the distance to the boundary of 
the marked vertex be less than D. In the regime z < 1/8, this distribution function is 
independent of z and identical to the two-point function of the Brownian map P,[7|,p!^ . 
This follows from the fact that, in the regime z < 1/8, the length of the boundary does 
not scale with n but remains finite at large n. This property can be measured as follows: 
as we already noticed, the generating function Wq counts quadrangulations with one 
marked edge along the boundary. To remove this marking, we must consider instead 
the generating function dz' {Wo{z') — l)/(2^') , so that the average half-perimeter 
reads 

Jo ? Is" 



From (|3.28| ), we immediately get that, at large n, 

so that {p)n{z) tends to the finite value 

, V . X z(zA(z))" Az , , . , . 

<P>"W-(iyF^- (l-8.)(l-VT^) 

Let us now consider the boundary-boundary distance statistics for ;2 < 1/8. When 
g 1/12, a sensible scaling limit is now obtained by keeping d finite in Ta- This is 
consistent with the fact that the perimeter itself remains finite and d < p obviously. For 
finite d, we now have the expansion 

d+1 {d+l){d + 2) 
= v^e + Oie) (4.26) 

while the expansion ( ^78| ) still holds. Using (|3.32| ), we get 

A(zf(2(d+2)-(d+l)A(z))\ , , , _x , , 
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In particular, decays exponentially with d as exp {—d/^{z)) with a correlation length 

log{A{z) - 1) 

If we now wish to compute the large n behavior of the term T^lgn , we have to extract 
the singular part of T^, which requires continuing the expansion ( ^4.27| ) up to order e"^/^. 
This yields a rather complicated and non-universal expression which exhibits the same 
exponential decay in the distance d. 

To conclude, the regime z < 1/8 is characterized by a perimeter which remains 
finite at large n and by boundary-boundary distances which also remain finite and are 
governed by a non-universal probability law with an exponential decay. This is to be 
contrasted with the bulk-boundary distances, which scale as n^/^ and are governed by 
the universal two-point function of the Brownian map. 

4-3. Scaling limit: the z > 1/8 regime 

Let us now discuss the situation z > 1/8, and more precisely l/8<2<l/4. In 

(2) 

this regime, we have a dominant singularity of the generating functions at Qcriti^) ^-nd 
we set 

g = gi^l{z){l~ue) (4.29) 
with e — s> 0. We then have the expansions 



. „,.(.) with .„,W^^^^-^-vW 



(4.30) 



2(1 -4z) 

,(2) 



Here a;crit(^) is the value of x, as given by (|3.4|), for g = g^Jti-^)- Keeping d finite, the 
generating function Wd has the expansion 

with Xcrit = a^crit(^) as above. We therefore find in this case a simple square root 
singularity in {g^}^ ~9)- This singularity (i.e the term proportional to -Jv e^l'^ in (^4.31| )) 
may be obtained alternatively from the behavior of VFd|zP at large p. Indeed, from (|3.25|) , 
we have at large p 



^ i-(/.+i-/.)E(2^+ix/^)'"\ 



k>l 



{4R)P (l-x^+^)(l-x'^+^)^ 

0Fp3/2 (l-xd+3)(l+xrf+l)2 
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(4.32) 



Upon summing over p with a weight z^^ the p-dependent prefactor gives rise when 

/2'\ 

ARz ^ 1 (which happens precisely when g approaches gcJti^) as in (|4.29| )), to a square 
root singularity — 2-\/l ~ 4Rz ~ —2 a/ (1 — 4^) / {8z — l)^/^ e^^'^, while x simply tends to 

its value Xcnt{z) at g = (/crin-^)- therefore recover the singular behavior in ( [4.31|) 
from the contribution of large values of p. 

As in the previous Section, we can consider quadrangulations with a fixed and large 
number n of inner faces. This is done again by performing a contour integral in g and 
setting 

9 = aSli^) {} + ■ (4-33) 

We may indeed write at large n 



-{.}^ ^^^^^^^-^^ / rfC^ e-M-} l + - • (4.34) 
2m g^+^ Tin J_^ 1 \ \n ' ' 



Using (|4.31| ) with e = 1/n and z/ = — C , we obtain that 



l^d o" ^ Q/o A/o 7 7: ^X777: ^XT77, ' ' 



0Fn3/2 V8z-l(l-xf+3)(l+xty)2 ' y^^'/' V8^-l- 

(4.35) 

As for log Wd , we have the expansion 



, T.. , /o(l-^crtt)(l+<rtt) \ /l-4^(l-Ctt)(l-<rtt) ^ 1/2 



with a singular part which can be alternatively read off the large p behavior 

(4i?)P (l-a;^+i)(l-a;'^+2) 



(4.36) 



(4.37) 



2v/^p3/2 (i + xrf+i)(i + x'^+2 
We immediately deduce the leading behaviors 



log Wd\ 



2v/^n3/2 V8z-l(l+a:t+')(l + 4t') 



logVF' -,(^SlW)-" /1-4. 



(4.38) 



20Fn3/2 V 8^ - 1 
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Fig. 8: Plots of the (non-universal) cumulative distribution function 
(pz{d) for z approaching the critical value 1/8 from above, namely z = 0.13, 
0.126, 0.1255 and 0.1251 (dotted plots from left to right). Beside each plot, 
we display the corresponding (universal) limiting scaling form (solid line) 
of Eq. (4.40). 



Taking the ratio log / log VFdIg" , we deduce the large n asymptotic expression 
4>z{d) for the probability that the distance to the boundary of the marked vertex be 
smaller than or equal to d in the ensemble of pointed quadrangulations with a boundary, 
namely: 

6 (d) - (l-(^c.u(.))^+^)(l-(xc.t(^))'^+^) .4 

' (1 + (Xe„t(^))^+l)(l + (Xc.t(^))^+2) ' ^ ' 

with d finite and Xcrit(^) as in ( [4.30|) . This function is expected to be non-universal. 
However, when z approaches the critical value 1/8, we have 



(t){d) ~ tanh^ {d(3{z)) with (3{z) = 2^3^ z - - , (4.40) 

and we expect that, except for the precise value of /3, the above scaling form is universal. 
The function (f)z{d) is plotted against its scaling form ( |4.40| ) for z = 0.13, 0.126, 0.1255 
and 0.1251 in Fig. 8. 



From the exponential growth with n of Wq \ g-n. , and from the general formula (|4.23| ) 
we now have at large n 



{pUz)-nxgSUz)U9Sli^))-' 

= nx ^^—^ (1/8 <z< 1/4) . 
1 - 4^; 
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More precisely, the probability to have a prescribed value of p is proportional to 
Wolgr^zpZ^ / (2p). Using the explicit form (|3.28|) , and expanding it at large n with p oc n, 
we find that, asymptotically, this probability tends to a Gaussian distribution peaked 
at p = {p)n{z) as above, and with width ^/n (1 — 4:z)/\/4z. 

Let us now consider the boundary-boundary distance statistics in the regime z > 
1/8. From the expression ( |3.32[ ) and the expansion ( |4.3U[ ), we see that a sensible scaling 
limit is now obtained by taking d large as 

d = De"^/2 (4.42) 

with D finite. Using fi = gB? ~ (1 — 4;^) / (122;) at leading order in e and fd+i — > 3;crit(^) 
for large d, we obtain that, in the scaling limit ( [4.42| ), 



_lJe-"V— V--. (4.43) 

We can again consider the fixed n ensemble by performing a contour integral in g. 
Taking d — Dn^/"^, we now get 



,(2) 

Li 

\ Z I I 1 



\ z ) V 8^ — 1 



(4.44) 



To have a proper probability density, we must multiply this quantity by the infinitesimal 
step rt}l'^dD and normalize it by the generating function of quadrangulations with two 
marked edges on the boundary, given at large n by 




(4.45) 



2. -TVol." - 2. - I — ^=-3^^ ^— ^ 

v/^nV2 V 1 - 4z V 2z 
We obtain finally the probability density 

1 — A.Z n2 l-4z 

Pbound.(^) = 215 e-^ — . (4.46) 

oZ — i 

Here Pbound. {D)dD measures the probability that the two marked edges on the boundary 
be at a distance in the quadr angulation in the range [-D, d + dD], in the ensemble of 
quadrangulations with two marked edges on the boundary. It is natural to measure the 
boundary-boundary distance d in units of the square root of the average half-perimeter 
{p)n{z), as given by ( [1.41|) . This is done by introducing the variable 



n 



^ - ^= = dJIzII , (4.47) 



l-4z 
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which remains a finite quantity in the scahng hmit. The probabihty density for the 
variable S foUows simply from (|4.46| ) and reads 




2 



/5bound.(5) = \ — Pbound. Si — = 26 6 ' . (4.48) 



This probability density is independent of z and, as announced in Section 2, is identical 
to the two-point function of the Brownian Continuum Random Tree [ESI, i.e a simple 
Rayleigh law. 



To conclude, the regime z > 1/8 is characterized by a perimeter which is pro- 
portional to n at large n and governed by a Gaussian law peaked at its average value 
( [4.41|) . The distance of a point in the bulk to this boundary remains finite and gov- 
erned by the non- universal distribution (pz{d) above. The boundary-boundary distance 
is of order n^/^ and, when measured in natural units given by the square root of the 
average perimeter, is characterized by the universal two-point function of the Brownian 
Continuum Random Tree. 

4-4 ■ Scaling limit: the critical regime 

Let us finally discuss the vicinity of the transition point z — ^crit = 1/8. A sensible 
scaling limit is now obtained by setting 

- = ^(l-^Be^/^) (4-49) 

o 

where d may now stand for both the bulk-boundary distance (as in Wd) and the 
boundary-boundary distance (as in Td). Using the expansions ( ^?6|) and the expan- 
sion 

W = 2(^1- ^^B + ^e'/^ + --^ , (4.50) 

we get 

Wd = 2(1- H{D; //s) + ■ ■ ■) (4.51) 

where 

HiD; ^, ^b) = fiD; ^) + " x^/^ (4.53) 

with /(-D; n) given by ( [4.7|) . The scaling function 'H{D] fx, hb), obtained here by a direct 
scaling limit of the discrete expression ( |3.12| ) for Wd, can be obtained alternatively as 
the solution of a non-linear differential equation as follows: using the expansion for Rd 
(as obtained for instance via ( [4.13|) for p = 1 since Rd = Wd\z^) 



Rn 



2 (^1 - J^(L>; ;u) + . . .) (4.53) 
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with T{D; n) given by and expanding the recursion relation (|3.6| ) at order e^/^ 

with Wd as in (|4.51|) , we get the equation 

doniD; fi, fiB) - n''{D; ^, ^b) + J'iD- fi) + = . (4.54) 

The expression ( [4 .521) above for 7i(L>; /x, (xb) is then the unique solution of this equation 
satisfying H{oo] fx, hb) = \/Ti^B~+~\/^, as required by (|4.50|) . We may also relate our 



expression to the result of Refs. and 0, by introducing the quantity 

g*iD;n,f,B)=dDd^^niD;f,,f,B) • (4.55) 

This is indeed the continuous counterpart of the generating function considered in 
Refs. and (in the slightly different context of triangulations) corresponding in 
our language to pointed rooted maps with a boundary where the origin-boundary dis- 
tance has a fixed value De~^^'^ (hence the operator do) and where the root edge lies 
anywhere on the boundary (hence the operator 9^g). With our explicit expression for 
7i, it is easy to check that Q* satisfies: 

doQ* = -2d^, {ICg*) with /C = Kill, ^ib) = (^I^B - ^) ^^B + v^ . (4.56) 

This is precisely the equation used in Refs. and to determine Q* and the two-point 
function. 

We have finally the expansion 

log Wd = log(2) -n{D;f,,iiB)e^/^ + --- (4.57) 

The regime ( [4.49|) corresponds to typical values of the perimeter of order e~^/^. 
Rather than fixing hb, we may alternatively work directly with a fixed value of the 
half-perimeter p being of order e~^/^, i.e consider our fixed length generating functions 
of Section 3.3 and set 

p = Pe-i/2 (4.58) 

Upon setting k = K t~^^'^, the expression ( |3.25| ) translates into 
e^'^niD.P-.l^) with 

(4.59) 
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The scaling functions Ti. and Ti are then simply related by 

/ 1 f°° ( e~^^ - \ 

n{D;fi,fiB) = ^^^B + VJ^+^J^ dPe-''-^(-j=^-n{D,P;fi)j . (4.60) 

For small P, we have in particular 

n{D, p- ^,f^' ^=^(1 - v^P + ■ ■ ■) (1 + Pi^VJ^ - V\D- ^)) + ■ ■ ■) 

1 (4.61) 

Note that, as could be expected, this small P behavior matches precisely the large p 
behavior of VFdl^p obtained in the regime 2; < 1/8 for finite p, as given by ( |4.13|) . 
For large P, we get instead 



^^^tanh^ 



(4.62) 



Again, as expected, this expression matches precisely that obtained for Wd\zp in the 
regime z > 1/8, as given by ( |4.32| ), upon considering the scaling limit p — Pe~^/^, 
P = 2(1 - 7/Ie^/2), d = De-^/^ and x = 1 - V6/i^/^e^/^. As for log W^, since, up to a 
factor 1/2, it has the same singularity than Wd, we have 

,m^.n{D,P;^). (4.63) 



This behavior can alternatively be obtained by taking directly the scaling limit of (|3.26|) . 

We can now turn to the fixed n ensemble, with n large and in the critical scaling 
regime 

p = P (4.64) 
with P finite. We then get a cumulative distribution function 

/OO (; 
- e-^' n{D,P;-^^) , (4.65) 
-00 ^ 

which measures the probability that a vertex chosen uniformly at random in the quad- 
rangulation be at a rescaled distance less than D from the boundary. This distribution 
function is plotted in Fig 2 for P = 0.01, 0.1, 0.5, 1.0, 2.0 and 5.0. For fixed P and 
small D, we have the expansion 

^{D,P) = ^PD^ --{P^-1)D'^ + --- (4.66) 
4 8 
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Fig. 9: Plots of the cumulative distribution function P) as a func- 
tion of the scaling variable Dy/P in the regime of large values of P (thin 
solid lines) and their comparison with the limiting scaling form of Eq. (4.67) 
(thick blue line). The plots represented here are for P = 1., 2. and 5., from 
left to right. 



For fixed D and small P, using ( [4.61 ), we immediately see that ^{D,P) ^{D) 
as expected. This property is illustrated in Fig. 2. Note that the small D and small P 
limits do not commute. 

On the other hand, when P is large, using ( [4.62| ) and evaluating the integral over 
^ by a saddle point estimate, we deduce that 



$(D,P)''=.°°tanh2 ( Vp 



(4.67) 



This property is illustrated on Fig. 9 for P = 1., 2. and 5.. 

Concerning the boundary-boundary distance, using the expression ( |3.33| ), we get 
in the scaling regime (|4.58| ) the limiting expression for Td\zP'- 



-dlzP 



8 e-VJ^P-^ 

3^ ^pr/, {{2D^-3DP) + {4D^P-2P^)f{D;f,)+2DP^f{D;fi)} 



(4.68) 

It is natural to measure the boundary-boundary distances in units of y/p, i.e write 
d — or equivalently 

D = 6Vp . (4.69) 
In the variable S, Eq. ( [4.68|) translates into 

-^-ge |(2^3-3J) + (4J^-2)yP/(VP(5;/x)-f2(5 (v^/(v^5;^)) 

(4.70) 
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At small P, we have the expansion 

-^^rij^Y^^^^ — 10 — 105 ^ ^ +---| 

(4.71) 

while at large P, we have 

Zik-""4y^7.£-^{2.e-%...} (4.T2) 

We turn finally to fixed values of n and p, in the critical scaling regime n 00 
with the ratio P — p/'u}/'^ fixed. Using ([4.70D with e = 1/n and n = — multiplying 
by the elementary step n^^^dD and normalizing by 2pWo\gnzP, which behaves as 

we get a critical probability density for the rescaled boundary-boundary distance D — 



PhoundXD.P) = -ije-^'/^{(2D3-3I}P) + (4D2p-2p2)ai(D,P)+2DPV2(D,P)} 



3P4 

y2 , 



with a,{D,P) = / d^ i e-^'+'^^f{D;-e) 

VttP 1 



(4.74) 

As before, it is natural to measure the boundary distances in units of y^, i.e write 
d = 5y/p or equivalently D = 5\fP. We find for 5 a probability density 



Pbound.(5,P) = ^e-'" {(253-35) + (452-2)ai(5, P)+25a2(5, P)} 
with ai(5,P) = /"^ ^ e-^^+'^^/(5v^;-e^) 

^2(5, P) = ^ / di i e-^^+^^^/^(5v^;-a ■ 
J-00 1 



(4.75) 



This probability density is plotted in Fig. 3 for P = 0.5, 1.0, 1.5, 2.0, 3.0, 5.0 and 10.0. 
At small P, we find in particular from ( [4.71| ) that 



Pbound.(^,^)''-^°T^e-'5' (355 + 28^3 + 1255 + 3^7) _ (4 75) 

105 
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At large P, we have the simple result 

Pbound.('^,^)'=^^2 5e-'''' . (4.77) 

Repeating the above analysis for the refined generating function Td{s, s') of 
Eq. ( |3.34| ), we have access to the refined probability PhoundX^^'^^ P) ^^^^ ^^e vertex 
of the boundary at rescaled distance 2Pu {0 > u > 1) along the boundary from a given 
vertex (chosen uniformly at random on the boundary) be at rescaled distance 5\/P in 
the quadrangulation from this vertex. This probability density reads 



1 g 4ii,{l— u) r 

ho...XS,u,P) = _-5^,^_^{(,^-2„)(^^-2(l-„)) ^^^^^^ 

We have in particular, for small P, 

Pbound.(5,^.,P) ^°^ ^5/2(i_^)5/2 l4^{70 + 21^^+3^^-42Kl-z.)(J^ + 5)} (4.79) 
while, for large P, 

p 1 5 e 4ii(i-u) 

Pbound. (<S, U, P) „3/2(l_„)3/2 ' (^'S") 




Fig. 10: The average value {5{u)) (measured in units of \/P) for the dis- 
tance in the quadrangulation of two boundary points at rescaled distance 
u (measured in units of 2P) along the boundary, here for P = 0.5, 1.0, 1.5, 
2.0, 3.0, 5.0 and 10.0 (thin lines from top to bottom). The plots interpolate 
between two limiting laws (thick lines) : the non-trivial (but universal) law 
of Eq. (4.82) for P ^ and the simple semi-circle law of Eq. (4.83) for 
P ^ oo. 
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A simpler characterization of the distance in the quadrangulation of two boundary 
vertices at distance 2Pu along the boundary is through the corresponding average value 



oo 



{5{u))p = d5 5 p^oundiS, u, P) . (4.81) 







This average value is plotted in Fig. 10 for P = 0.5, 1.0, 1.5, 2.0, 3.0, 5.0 and 10.0. For 
P small, we get in particular 



{S{u))p''^'i^J^^^^—^\35 + 21u{l -u)+ 36u'{l - u)^} (4.82) 

105 V TT I J 

while, for P large, we have 



(^(„)),-=.»4^^ . (4.83) 



To conclude, quadrangulations with a boundary display an interesting scaling be- 
havior when both the number n of inner faces and the length 2p of the boundary become 
large, keeping the ratio P — p/n^^"^ fixed. In this regime, the bulk-boundary distances 
scale as n^/'^ and are characterized by a universal distribution function $(-D, P) which 
interpolates between the (cumulative) two-point function of $(-D) the Brownian map 
for small P to a simple tanh'^ function in the variable D\fP at large P. As for the 
boundary-boundary distances, they scale as the square root of the perimeter and are 
characterized by a universal probability density /5bound.('^5 which interpolates between 
the two-point function of the Brownian Continuum Random Tree when P is large to the 
small P expression ( f4.76| ). Note that this latter formula, even if it looks rather involved, 
is expected to be universal, and so are the formulas (|4.79|) and ([4.82|). 



5. Self-avoiding boundary 

5.1. Generating functions 

So far, we considered quadrangulations whose boundary may contain separating 
vertices or edges, corresponding to vertices or edges encountered several times along 
the contour. We may instead consider quadrangulations with a self-avoiding boundary, 
i.e demand that the 2p vertices (and consequently the 2p edges) along the contour 
be all distinct. We shall denote by Wd the generating function for quadrangulations 
with a self-avoiding boundary having a marked vertex at distance smaller than or equal 
to d from the boundary, and as before with a marked "closest edge", i.e a boundary 
edge incident to a vertex at minimal distance from the marked vertex and oriented 
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Fig. 11: A schematic picture of the Eq. (5.1). In any rooted quadrangula- 
tion with a boundary, the root edge selects a particular irreducible compo- 
nent (in magenta) which may be either a single edge (a) or an irreducible 
component with non-zero area and self-avoiding boundary of perimeter 
2pQ. The other irreducible components may be reassembled into (possi- 
bly empty) quadrangulations with a generic boundary attached to each of 
the two endpoints of the edge (case (a)) or to each of the 2po boundary 
vertices of the selected irreducible component (case (b)). These attached 
quadrangulations are naturally rooted and each of them gives rise to a 
factor Wo{g, z). 

counterclockwise around the bulk of the quadrangulation. As before, when (i = 0, these 
markings reduce to the choice of a boundary edge oriented clockwise around the bulk. 

In Wdi we weight each configuration by a factor g per inner face and a factor \/~Z 
per edge of the boundary. For convenience, we decide to also add to the configurations 
counted by Wa the empty configuration (weight 1) and the configuration with p = 1 
and 7T, = (weight Z) corresponding to a single edge embedded in the external face, 
although the boundary is not self-avoiding in this case. In other words, we take the 
convention that VFdIgO = 1 + Z in the following. With this convention, it is easy to 
check that we have the combinatorial identity 

Woig,z) = Woi9,Z) with Z = zW^{g,z). (5.1) 

Indeed, considering the root edge in any (non-empty) configuration counted by Wq, with 
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the external face on its right, it belongs to some irreducible component which is either a 
quadrangulation of non-zero area, with a self-avoiding boundary of length, say 2pQ with 
Po > 1, or a single edge (po = 1) if the root edge is a separating edge (see Fig. 11 for 
an illustration). The configurations in the latter case are counted by zWQ{g^z), while 
those in the former case are counted by (Wq — 1 — Z)\zpo x [zWq {g,z)Y° since we 
may attach to each vertex of the self-avoiding boundary of the irreducible component a 
configuration counted by Wo{g, z). Summing over po, we end up with the relation (|5.1| ). 




Wo(g,z) 



Wo(g,z) 



Fig. 12: A schematic picture of the Eq. (5.2). In any pointed-rooted quad- 
rangulation with a boundary counted by Wd — Wq, the origin vertex lies 
strictly in the bulk of some particular irreducible component (in magenta), 
with self-avoiding boundary of perimeter 2po. The other irreducible com- 
ponents may be reassembled into 2pQ naturally rooted quadrangulations 
with a generic boundary attached to the 2po boundary vertices of the se- 
lected irreducible component. One of these quadrangulations contains the 
originally marked closest edge, inducing a natural splitting in two parts (in 
dark red) which results in an extra factor Wo{g, z). 

For (i > 0, we have the slightly different relation 

Wd{g,z)-Wo{g,z) = Wo{g,z)(Wd{g,Z)-Wo{g,Z)') with Z = zW^{g,z). (5.2) 

Indeed, Wd — Wq counts configurations with a marked vertex strictly in the bulk of the 
quadrangulation, at distance less than or equal to d and with a marked closest edge. 
The marked vertex belongs to one particular irreducible component and its distance 
to the whole boundary is equal to its distance to the self-avoiding boundary of this 
component. The most general configuration is then obtained by attaching to each vertex 
of the self-avoiding boundary of the irreducible component a configuration counted by 
Wo{g, z), leading again to the effective weight Z = zWq in the generating function for 
the irreducible component (see Fig. 12 for an illustration). Finally, the marked closest 
edge on the original configuration starts form a particular closest vertex on the self- 
avoiding boundary which in turns selects a closest edge on the self-avoiding boundary, 
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leading eventually to Wd — Wq- The marking of the original closest edge induces a 
splitting in two parts of the contour of the configuration attached to the marked closest 
vertex on the self-avoiding boundary, resulting in an extra multiplicative factor Wo{g, z) 



Fig. 13: A schematic picture of the decomposition of a pointed quad- 
rangulation with a self-avoiding boundary, whose origin lies strictly in the 
bulk and with k closest vertices on the boundary (here k = 7), into k 
slices, as defined in the text. From each closest vertex on the boundary 
(red dots), we draw the leftmost geodesic path to the origin (blue circle). 
Note that these geodesies cannot cross but may merge before reaching the 
origin (as in V3). This splits the quadrangulation into k triangular slices 
whose depths are generally different, with the maximal depth being equal 
to the origin-boundary distance. 

It is more natural to consider configurations without a marked closest edge. In 
the case of a non self-avoiding boundary, this was done by changing Wd into logWd- 
Here, we find that suppressing the marking of the closest edge results in replacing the 
generating function Wd by a modified generating function 



A proof of this relation is sketched as follows. It is equivalent to exhibit a generating 



in (13). 






(5.3) 
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function Sd such that: 

Wd - (Wo - 1) = 

fc>0 

- . (5.4) 

k>i ^ 

CombinatoriaUy, both left hand sides are generating functions for pointed quadrangu- 
lations with a self-avoiding boundary, such that the origin-loop distance is between 1 
and d, with an additional marking of a closest boundary vertex in the first case (note 
that the marking of an edge or of a vertex are equivalent for self-avoiding boundaries). 
Clearly, we shall look for a decomposition according to the number k of such closest 
vertices. A possible way is to decompose the map into k "slices" by cutting along each 
leftmost geodesic from a closest vertex to the origin, see Fig. 13. Each slice is itself a 
quadrangulation with a self-avoiding boundary, and has three distinguished vertices on 
the boundary: two of them, say Vi and V2, correspond to two clockwise consecutive 
closest boundary vertices in the original map and the third, say vs, is the point where 
their outgoing leftmost geodesies merge (note that ^3 is not necessarily the origin). 
Furthermore, the slice contour consists of three segments {V1V2), {v2Vs) and {vsVi) of 
non-zero length, and the following properties hold within the slice: 

- {vsvi) and {v2Vs) are geodesies of the same length £, hereafter called the depth of 
the slice, with 1 < i < d, 

- there is no other geodesic from vi to f 3 (due to our convention of cutting along 
leftmost geodesies), 

- all vertices on (f 1^2) distinct from vi and V2 are at distance strictly larger than £ 
from V3. 

Let Sd be the generating function for slices satisfying the above characterization, with a 
weight g per face and ^/Z per edge of (^1^2)- Then, we have a bijective decomposition 
of the original map into a linear or cyclic sequence of slices, depending whether a closest 
boundary vertex is marked or not. Note that the slices obtained in the decomposition do 
not necessarily have the same depth, and that the maximal depth is equal to the origin- 
loop distance. Omitting further details of the proof, the identities ( |5.4| ) are established. 

The generating function for pointed quadrangulations with a self-avoiding boundary 
with origin-boundary distance equal to d is 

(5^ = VT^-VTd-i, Go = K = Wo-l , (5.5) 

which yields the formula ( |2.12| ) announced in Section 2. 

Let us now derive explicit expressions for Wd and W^. From (|5.1|) and from the 
general expression ( p.lSj ) for Wq, we immediately deduce that 

Wo{g.Z) = W{l-h (W-l)) , (5.6) 

where 

W = W{g, Z) = W{g, z{g, Z)) with z{g, Z) given implicitly by 

Z = z{g, Z)W^{g, z{g, Z)) = z{g, Z)W^{g, Z) ^^'^^ 
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Writing these equations as 

W = l + z{g,Z)RW^ 



Z = z{g,Z)(wil-hiW-l))) , 



and eliminating z{g, Z), we deduce the relation 

„ W-1 



R 



l-f,{W-l)r (5.9) 



which determines as a function of g and Z. As for Wd for (i > 0, we deduce from 
( ^■2|) and from the general expression (|3.12| ) for Wd that 

Z) = 1 X ^-iW-l)fd+i ^ 

"^'^ l-f,{W-l) 1-{W-I)fd 



with W related to Z as in ( p.9| ) above. Note that the quantity Wd depends on Z only 
via W but that, unlike generic boundaries, W is not the limit of Wd for (i ^ oo. Finally, 
we have 

W:,{g, Z) = log ( 1 X 1 - (W- l)fd+i 1 + _ 1 . (5.11) 

As before, Wd and W^ are Lagrangean generating functions and we may extract 
an explicit expression for their Z^ term. Writing ( p.9|) as 

Z = ^{l-fiwr (5.12) 

upon introducing the quantity 

w = W-l, (5.13) 

we may easily transform any contour integral in the variable Z into a contour integral 
in the variable w, namely 

' n'i'i^Ar..^^M^i\- (5.14) 



2i7r ZP+i I J / 2i7r wp+^ (1 - fiw)^P+^ 

Using (|5.6|) , which we write as 

Wo = il + w)il-fiw) , (5.15) 

and /i = g , we obtain 
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for p > 1, and Wo\zo = 1. Extracting the g'^ term of the right hand side, we deduce the 
equivalent formula 

W\ ----- ^^P^- C^n + P -1)1 

valid for all n > 0, 1 < p < n + 1. 

As for Wd, we get upon expanding ( |5.1C1| ) in w = — 1 the expression 

Wd = ^— ^ I 1 - * ifd+i - fd) $^(* /d)'-' \+Wo-l (5.18) 
from which we deduce 

w,w, = (sii')^ -2(^3_f J (S^^^p^/,)-. 

- r.-,p3^J (3p)! _ . . (3p-fc)! (2/c + l) 

" ^ |p!(2p+l)! ^^^+' ^^^^^ (p-A))!(2p+l)! {gR'^)'^ ^^"^ 

+ VFolzp - 5p,o ■ 

(5.19) 

By a similar argument, we have 

-^1- - <««^'' { W - ' t - } (5.20, 

+ WqIzp 

for p > 1. 

5.^. Critical behavior and scaling limits 

As before, the generating functions Wd and have a first singularity at g = 1/12, 
irrespectively of the value of Z. A second singularity comes from W which, from (|5.9| ), 
is singular when 

27 

-gR^Z = l. (5.21) 

This equality may occur only when Z > 2/9 and it defines a critical line g = gcrit{Z), 
with 

gcritiZ) = ((4 + 9Z) VieZ + 9Z2 - 9Z(4 + 3Z)) (Z > 2/9) . (5.22) 
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For Z > 2/9, we have ^crit(-^) < 1/12 so that this new value determines the radius of 
convergence in g of Wd- As before, we have a change of determination of the radius of 
convergence at the critical value 

2 

ZcTit = - 5 (5.23) 

corresponding to a transition between two differently behaved regimes. 
To study the regime Z < 2/9, we set as before 



and we use the expansion (|5.19| ) to get 



(6.24) 



- (2f(D-, e'/^ ± 3'-^-i|p^(2* + 1) } + • • • (5.25) 



Here we have used the identity 



y 3^^-%^(2^ + 1) = 3 , I'f^' „ (5.26) 

^ (p-/c)!(2p+l)!^ ' (p-l)!(2p+l)! ^ ' 

obtained by writing the summand in the left hand side as ak — ctk+i with ak = 
^k^3p-k-iy gijj^iiarly, we have 



W, 



d\ZP 



for p > 1. 

Finally, from ( |5.16| ), we have: 

'^"l-^KD'il^fT'^+^W' (5.28) 

for p > and Wo\zo = 1. Summing (|5.25|) over p with a weight Z^, we obtain 
Wd = Ao{Z) + A{Z) - 3ZA'{Z) T{D- ^) e^/^ + • ■ ■ 

2 11 ni 9Z 



with Ao{Z)=2y(-] ZP^^^^ = -(-!+, F, 

v^y pK^P-iy. 3 V 



3' 3 



1 



i(Z)=y f-V Z^-JML_ = ^/_sin I -arcsin | 3i 
p!(2p+l)- 



L)! 



(5.29) 
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while 

W', = Ao{Z) + log A{Z) - T{D- ^) + 

Going to a fixed n ensemble, we deduce from this expansion that 

12'" ^ZA\Z) 



(5.30) 



dig' 



7rn3/2 A{Z) 



/oo 
-oo 



W'Lr^ - ^ ^ where W' = lim 



(5.31) 



at large n so that, taking the ratio of these quantities, we get the same asymptotic 
distribution function ^{D) as in Section 3 for the rescaled distance D. In the Z < 2/9 
regime, the length of the boundary does not scale with n and quadrangulations with a 
self-avoiding boundary stay, in the scaling limit, in the universality class of the Brownian 
map. 

To study the regime Z > 2/9, we must now set 

Z = g^,it{Z){l-i>e) (5.32) 
and keep d finite. We introduce the notations 

Xcrit(^) = ^(27Z-8 + 9V^(16 + 9Z) 

-V6^243Z2 + 144Z-32+3(27Z-8)V"Z(16+9Z)) (5.33) 



^crit _ ~ fy\ I - (XcTitiZ))'^ 



1 - {Xcrit{Z)) 



d+2 • 



Here 5crit(-^) is simply the value of x when g = gcvitiZ)- With these notations, we have 
the expansion 

W= fl + ^l - (7(Z)v^6V2 (5.34) 



where C{Z) is some function of Z (with no singularity for Z > 2/9) which we do not 
make explicit. We then find the expansion 



^ _ 4+21/r^ ^ 3 

Igjcrit 2 (3/crit _ Jcrit) 

c>m./~,i/2f 4+i5/r^ 9ff' ifY' - m\mY' - ft') 

y ) S 12 4 (3/crit _ Jcrit-)2 
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(5.35) 



40 60 80 100 



Fig. 14: Plots of the (non-universal) cumulative distribution function 
(pzid) for Z approaching the critical value 2/9 from above, namely Z = 
0.23, 0.225, 0.223 and 0.2227 (dotted plots from left to right). Beside each 
plot, we display the corresponding (universal) limiting scaling form (solid 
line) of Eq. (5.38). 



and the similar expansion 



crit 



9/r 



crit 



^ , ,3 (3/r^ 

+ log 7^ 



i"crit 
•I ( 



2 crit ycrit^i 



-(7(Z)v^ei/2 



2 + 3/f 
6 



crit 



+9(/r')' 



/ /crit 



i-crit 



(5.36) 



(3/f 



crit 



/crit 
J d+1' 



In the fixed n ensemble, we deduce that the ratio W^\g^ /W'^\g^ tends at large n to the 
cumulative distribution function for d: 



2 + 3/f"' (3/i"" - /5"')(3/f"' - /J"' 



' d+l! 

2\ ( ( 2+x ^d\( 2+x 



27x(l + x + x2) 



^ (2 + x)2(i + 2x)2|^ (r^ + ^'^)(r^ + ^'^+^)J ^^'^ ^ - • 

(5.37) 

Note that this distribution is different from the distribution ^^((i) of Eq. (|4.39|) obtained 
for a non self-avoiding boundary. When Z approaches the critical value 2/9 however, 
we have 

3 r, 2 



\z{d) ~ tanh^ (rf/3(^)) with f3{Z) = - J Z - - , (5.38) 
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and we recover the universal scaling form of Eq. ( [4.40|) . The function (f)z{d) is plotted 
against its scahng form (|5]3|) for Z = 0.23, 0.225, 0.223 and 0.2227 in Fig. 14. 
Finally, to study the regime Z ~ 2/9, we can set 



Z = lil-~,se^/') (5-39) 



2 

if we work with a fixed value of Z, or alternatively replace the second equation by 
if we work in a fixed p ensemble. In the first ensemble, we have the expansion 



Wd - (VTo - 1) = ^ (l - e'/'^ {D; ^, Ab/3) + • ■ ■) 
VT^ - {Wo - 1) = log - e'/^n (D; ^, /is/S) + ■ 



(5.41) 



with Ti. as in ( ^ ), while Wq = 4/3+0(ei/2). We thus recover the same scaling function 
as that obtained for non self-avoiding boundaries, with a simple renormalization of the 
"boundary cosmological constant" by 1/3. 

In the fixed length ensemble, we obtain from ( p.l9| ) the scaling behavior 



(9/2)P 4 

with Ti. as in (|4.59| ). Except for the trivial factor 9/4, we recover again the same scaling 
function as that obtained for non self-avoiding boundaries, with now a renormalization 
of the perimeter by a factor of 3. As for W^, since it has the same singular behavior as 
Wd up to a factor 3/2, we deduce 

e3/4 ^n{D;3P;^) . (5.43) 
(9/2)P 2 ^ ' '^^ ^ ^ 



This behavior can also be obtained directly by taking the scaling limit of (|5.20|) . 



After going to the fixed n ensemble and taking the appropriate ratio W^\g^zp /W^\gr^zp , 
we now find that 

$ .clf-avoiding (D, P) = $(D, 3P) (5.44) 

boundary 

for the distribution function for D in the fixed p ensemble. 
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Note that the ratio W^lzp /W^\zp tends to 3/2 which may be interpreted as the 
asymptotic average number of boundary vertices closest to the origin in the case of a 
self-avoiding boundary. This number is independent of D in the scaling regime. 

To conclude, large quadrangulations with a self-avoiding boundary behave essen- 
tially as large quadrangulations with a non-self-avoiding boundary, as far as the bulk- 
boundary distance is concerned. For boundary lengths of the order n^^"^, quadrangu- 
lations with a self-avoiding boundary of length P behave as quadrangulations with a 
non-self-avoiding boundary of length 3P. This suggests that, in quadrangulations with 
a non-self-avoiding boundary, only one of the irreducible components is macroscopic 
and has a boundary of length equal to 1 /3 of the total perimeter, while all the other 
components are microscopic but altogether, the lengths of their boundaries represent 
2/3 of the total perimeter. 

6. Self-avoiding loops 

We end this paper by a study of the statistics of distances in quadrangulations with 
a self-avoiding loop. As already discussed in Section 2, a self-avoiding loop is a closed 
path made of consecutive edges of the quadrangulation, which is simple, i.e visits any 
vertex at most once. For convenience, we shall suppose that the loop is oriented. Upon 
cutting along the loop, we obtain two quadrangulations with a self-avoiding boinidary, 
constrained to have the same perimeter. The orientation allows to distinguish this two 
pieces as left and right. We can now express a number of generating functions for this 
problem in terms of the generating functions found in Section 4. 

A first simple generating function is that of quadrangulations with a self-avoiding 
loop with a marked vertex on the loop. It reads immediately 

p>i 

where g is the weight per face while ^yy is a weight per edge of the loop (the length of 
the loop is necessarily even). 

Another interesting quantity involving the distance is the generating function for 
quadrangulations with a self-avoiding loop with a marked vertex at distance less than 
or equal to d, lying to the right of the loop (or possibly on the loop itself) . It reads 

^d{9,y) = ^yWo(^,Z)|z.W^^(^,Z)|z. (6.2) 
p>i 

if the configurations are counted with their usual inverse symmetry factor. This 
can be seen by first noting that the generating function for the same objects 
with an additional marking of a closest vertex on the loop is obviously given by 
'^p^iy^Wo{g, Z)\zpWd{g, Z)\zp upon gluing the two boundaries in such a way that 
their marked vertices coincide. Removing the marked closest vertex, we see that a given 
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configuration is over-counted by a factor of k equal to the number of closest vertices 
on the loop. Dividing by this factor k amounts to replacing Wd by as in Section 4. 
Note that, clearly, the quantity is related to the quantity of Section 2 by 

Td = ^d-^d-i, To = Oq. (6.3) 

From the exponential growth Wq ~ (9/2)^ and ~ (9/2)^, we immediately 
deduce that a transition occurs now aX y = ycrit with 

2\^ 4 



For y < 4/81, the typical values of p contributing to fld{g, y) remain finite when g tends 
to its critical value 1/12. Setting 

q = — (1 — u e) 

12^ ^ ^ (6.5) 

d = De-'/^ , 
we may use ( p.27| ) and ( p.28| ) to write 
o| ./4y/ (3p-3)! V|0^4y (3p-1)! (3p-3)! r . .1/2, \ 

(6.6) 

for p > 1. Summing over p with a weight y^, we obtain 
Od = a(|/)-6(|/)^(I);;u) e^s + . . . 

4y ^ (3p-3)! / (3p-3)! , (3p-l)! 



with «(^)-2|:(-J ^^^1^(2-^^+ ^ 

^ V9y' ^ (p-l)!(2p-l)! pl{2py. ' 

which are functions of y with no singularity for y < 4/81. Going to a fixed n ensemble 
and considering large values of n, we immediately deduce from the singular behavior 
( |6.7| ) that the cumulative distribution function for the rescaled distance D = djirt}!^ 
from the marked vertex in the bulk to the self-avoiding loop is again, in the regime 
y < ycritj equal to the universal two-point function $(-D) of the Brownian map, as given 
by Eq. ([4.22| ). In this regime, the size of the loop does not scale with n and becomes 
negligible in the large n limit. 

In the regime y > 4/81, we expect that the dominant singularity corresponds to 
large values of p. In this case, we may use ( |5.16|) and ( 5.20 ) to write the large p behavior 




3/1 

(l-3/i)(2+3/i 




2 + 3/1 i3h-fd)i3h-fd+i: 

(6.^ 
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The line y{27gR^ / A)"^ = 1 defines the critical value gcritijj) of g for y > j/crit, namely 

9cvit{y) = gcvitiVy) (6.9) 

with ^crit(-^) as in Eq. ( |5.22| ). Upon summing over p with a weight y^, the p-dependent 
prefactor gives rise when {{27/4)gR^)^y 1 to a singularity 



(6.10) 

This singularity translates into an asymptotic behavior of the form {gcritiv)) jv? for 
O^lg" at large n, with a multiplicative factor proportional to the d-dependent coefficient 
in (|6.8|), taken at x = x(5'crit(l/)) = 35(v^)- Note that this d-dependent coefficient comes 
from W'^ only and we therefore recover in the fixed n ensemble (and in the regime 
y > 4/81) the same form (|5.37|) for the cumulative distribution function 4>z{d) for d 



provided we identify Z = ^ in this formula. 

The most interesting situation is when y is in the vicinity of 4/81. Alternatively, we 
may work in the fixed length ensemble and study the scaling limit ( |6.5| ) with moreover 



P 



. (6.11) 



In this limit, we may use the relation (|5.43| ) and the relation 



(9/2)P 0F(3P)V2^ ^ v^; y ) 



inherited from (|5.16| ) to deduce that 



Q I p-^Vf^P 

3^-— ^(1+3PV^)7Y(D;3P;//) 



r\-/ 



(81/4)P v^(3P)V2 

^e^n{D,P]fi) with 

(6.13) 

involving a new universal scaling function 'H{D, P; n). 



For small P, we have 



n{D, P; ^,f^' (1 - {3P)T{D; ^) + ■ • ■) (6.14) 



which matches the large p behavior of fld\yp in the regime y < 4/81. For large P, we 
have instead 



?i(AP;^)^~°°3-^^v^tanh2(^^^//4l)l . (6.15) 
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If we now turn to the fixed n ensemble, in the Hmit of large n, with the scaling 

p = P ^ (^Q jQ) 

we get a cumulative distribution function 



(6.17) 



which measures the probability that a vertex chosen uniformly at random in the quad- 
rangulation with a self-avoiding loop of rescaled length 2P be at a rescaled distance less 
than D from this loop. This is a new universal distribution which cannot be reduced to 
the distribution ^{D, P) of previous sections by a simple rescaling of P. It is plotted 
in Fig 4 for P = 0.01, 0.1, 0.2, 0.5 and 1.0. For small D and fixed P, we have the 
expansion 

, 9 , 3 1 + 9P2 + 162^-* , 
^(D,P) = -PD^ D^ + --- (6.18 

From ( |6.14| ), we see that ^{D,P) $(-D) at small P for fixed as expected. 
From ( |6.15| ), and via a saddle point estimate of the integral over ^, we find that, for 
P oo, the typical value of D is of order with the scaling form 

^{D,P) ~ tanh^ (^\/f^ ' 
i.e we recover the general form found in previous Sections. 



7. Discussion and conclusion 

To briefiy summarize our results, we have been able to find discrete exact expres- 
sions for the bulk-boundary and boundary-boundary correlators for quadrangulations 
with a generic boundary. In the case of a self-avoiding boundary, we have been only 
able to express the bulk-boundary correlator, leading eventually to results for a model 
of self-avoiding loop. From our discrete expressions, we identified three scaling regimes 
for which we gave the asymptotic behaviors. 



Most of the above results follow from the discovery of the "master formula" (|2.5|) for 
Wd, which is remarkable in itself. This is yet another manifestation of the still mysterious 
integrability of the equations governing distance statistics in maps, as already observed 
for the two-point and three-point functions. Here two levels of integrability are involved, 
first for the equation determining Rd, involving the parameter g only, and then for the 
equation ( |3.6|) determining Wd, involving z and Rd as an "external potential". In this 
respect, it is worth mentioning that we have a continued fraction expansion 

Wd = ^ (7.1) 

^ _ zRd+2 
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and furthermore, we have the conserved quantity 

Wd-gzRdRd+iRd+2WdWd+iWd+2 = W-gzR^W'' for aU . (7.2) 

This identity may be checked directly from the exphcit expression ( |3.12| ) for Wd- Alter- 
natively, using ( pTsP to expand both sides of this equation in powers of z, it is equivalent 
to an infinite number of conserved quantities involving Rd and g only, and it is easy to 



check that those correspond precisely to the conserved quantities found in Ref. ||28|| . 

Another remarkable fact is the strong similarity, already apparent at the discrete 
level, between the bulk-boundary correlators for generic and self-avoiding boundaries. 
Unfortunately, we have not been able to exhibit the same phenomenon for boundary- 
boundary correlators. Furthermore, had we gone over this problem, we are still far from 
understanding the distance statistics between two points lying on a self-avoiding loop: 
the loop and geodesies may cross each other, which prevents the decoupling of both 
sides observed in ( |6.2| ) for the bulk-loop correlator. 

Our universal expressions are properties of what could be called the Brownian 
map with a boundary. In this respect, we may wonder whether these results may 
be re-obtained in a purely continuous formalism, which would likely require a proper 
probabilistic definition of this object. At a more physical level, we notice that the 



expression ( 4.78 ) may be re-derived from the solution V{D,U;fj,) of some diffusion 



equation in a potential J^{D; n) as in Eq. ([4.11| ): 

Here V describes the law of the position D at time U of a particle diffusing in the 
potential in one spatial dimension. Discarding normalization factors, it reads: 

p-V^U 2 

V{D, U- fi) cx ^T^e- w {D^ -2U + 2UDf{D; ^)) , (7.4) 

with f{D;^) as in (|4.7| ). The expression ( [4.78| ) can be obtained by considering the 
quantity V{D, U ; fi)V{D, P — U; fx), setting // = — and performing the usual appro- 
priate integral over ^ to go to a fixed area ensemble, and finally introducing the rescaled 
variables d = Dj^fP and u = U/P. Heuristically, V{D,U; n)V{D, P - U; n) is the 
continuous counterpart of the generating function for Dyck paths of fixed length, con- 
strained to reach some prescribed height at some given step. The potential JF(-D; fx) is 
the continuous counterpart of the weight Rd attached to each descent d ^ d — 1 of the 
Dyck path. 

To conclude, let us list a few possible generalizations of our results. We expect 
integrability to survive when considering maps with faces of degrees other than 4, as 



found for the two-point function in We may also consider maps of higher genus 
and/or with several boundaries. For instance, Refs. ^ and give continuous results 
for surfaces with two boundaries at a prescribed mutual distance, for which discrete 
formulas may as well exist. Finally, introducing multiple boundaries may pave the 
way towards understanding the distance statistics in the general 0{J\f) loop model on 
dynamical random lattices. 
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